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ABSTRACT: Te02 crystals are used as bolometers in experiments searching for Double Beta Decay 
without emission of neutrinos. One of the most important issues in this extremely delicate kind 
of experiments is the characterization of the background. The knowledge of the response to a 
particles in the energy range where the signal is expected is therefore a must. In this paper we 
report the results on the response function of a Te02 bolometer to a's emitted by 147 Sm dissolved 
in the crystal at the growth phase. A Quenching Factor of (1.0076 ±0.0005) is found, independent 
of the temperature in the investigated range. The energy resolution on a peaks shows a standard 
calorimeter energy dependence: a [keV] = (0.56 ± 0.02) (0.010 ± 0.002) y/E[keV] . Signal pulse 
shape shows no difference between a and j3/y particles. 
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1. Introduction 

Bolometric detectors |||] are used in particle physics experiments to search for rare events like 
Neutrinoless Double Beta Decay and Dark Matter interactions. They are sensitive calorimeters 
operated at ~10 mK that measure the temperature rise produced by the energy deposited in particle 
interactions. 

An array of bolometers made of Te02 crystals has been used in the CUORICINO experi- 
ment [Q, to search for the Neutrinoless Double Beta Decay of 130 Te. A precise knowledge of 
the response function of these bolometers to each species of particles (j3/y and a) is important for 
the background rejection. Indeed, a sizeable fraction of the background in the Double Beta Decay 
Region of Interest (Rol), that for 130 Te is around 2.5 MeV [Q, is due to degraded a particles 
that lose part of their energy in the detector support structure and the rest in a single bolometer 
mimicking a signal event. 

While the response of Te02 crystals to y interactions up to 2.6 MeV is well known from the 
routine calibrations performed with Th radioactive sources, there are still relevant items that should 
be addressed for a particles. The most important are the Quenching Factor (QF), the relative 
energy resolution and, finally, possible signal shape differences with respect to 7/ /3 interactions 
particularly in the region where the signal is expected. 

The QF is defined as the amplitude ratio between the signal produced by an a particle and the 
one produced by an electron depositing the same energy in the detector. It is expected to be very 
close to unity in thermal detectors since any kind of energy deposition should be converted into 
heat. 

The a's that can be normally measured in the bolometer are either coming from decays inside 
the crystals or by interactions induced by external sources. Those coming from internal contami- 
nations (usually from U and Th chains) have energies higher than 4 MeV, well above the Double 
Beta Decay region of interest. 
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The QF for a particles in the energy region [5.7-8.8] MeV was measured in Te02 detectors [Q] 
using a 228 Ra a radioactive source. It was found to be 1.020 ± 0.005 (stat) ± 0.005(yy.tf ). No 
deviation from constancy of QF as a function of energy has been observed. Nevertheless this 
measurement is based on an extrapolation of the calibration function at energies well above the 
region attainable with a Th source (2.6 MeV), where non-linearities of the Te02 bolometers play 
an important role [Sj. 

In this work we studied the response of Te02 bolometer to an internal source producing 
monochromatic a particles with energy close to the Q- value of 13() Te Neutrinoless Double Beta 
Decay. We achieved this goal dissolving into a Te02 crystal a small amount of natural samarium. 
The naturally occurring isotope 147 Sm has an isotopic abundance of (15.0±0.2)% [JTTJ] and an half 
life of 1.06- 10 n y [11 J. To have areasonable a decay rate (~count/hour) the crystal should contain 
only a few micrograms of natural Sm. 

The a particles are emitted in the 147 Sm — s- 143 Nd +a transition. The Q-value of this reaction 
is shared between the emitted a particle and the recoiling nucleus. Since the decay is contained 
within the crystal, we measure the full transition energy of (2310.5 ±1.1) keV [11, O]. 



Beside the obvious advantage of measuring the a response in the actual Rol, an additional yet 
extremely relevant feature is present and could be fruitfully exploited. Bolometers, in fact, require 
a continuous monitoring of the operating temperature. Even tiny drifts can induce variations of the 
thermal gain, thus spoiling the energy resolution. To correct for this effect, a Si resistor (heater) 
glued to the detector surface produces heat pulses by Joule dissipation which are very similar to 
particle induced pulses and the temperature drift can be corrected offline on the basis of the mea- 
sured heater pulse amplitude variation []|, [l(J. Heaters however can experience electronic failures 
or grounding problems, resulting in changes of their pulse amplitude not related to temperature 
drifts. This calls for a frequent re-calibration of the entire system that is obviously a long time 
subtracted to physics data taking. 

The presence of a long-lived monochromatic a line in the proximity of the Rol could take the 
role of the heater and allows less frequent calibrations. 

The outline of this work is the following: in Section ^| general properties of Te02 crystals and 
the doping process with natural Sm are described, in Section || details of the experimental setup are 
given while results are summarized in Section |j. 

2. Doping of Te02 crystal with an a-particle emitter 

Two methods are currently used for growing Te02 crystals: Czochralski and Bridgman. The growth 
is rather difficult, special temperature gradient conditions as well as pulling and/or rotating rates 
being needed in order to obtain high quality crystals. Further post-growth thermal treatments are 
applied to the as grown ingots aimed at quenching the tendency of Te02 single crystals to cracking 
caused by high anisotropy of thermal coefficients. Inhomogeneities in the crystal due to the incor- 



poration of Pt dissolved from the crucible is another problem frequently reported | |14| ] and careful 
control of convection currents in the melt is needed to avoid the incorporation of the gas bubbles to 
the crystals, especially in the case of Czochralski growth. 

Te02 crystals do not easily permit impurities inclusion in the lattice. The paratellurite struc- 
ture, especially the asymmetric covalent Te-0 bonds, limits the incorporation of the dopants to 
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extremely low levels which results in very low segregation coefficients (k~10~ 2 ). To avoid pre- 
cipitation or aggregation of foreign ions during crystal growth only small amounts of dopant can 
be added to the melt and the resulting built-in concentration is often below the limit for chemical 



analysis [14]. Only a few successful doping of Te02 are reported in literature for ions like Fe or 



Cr [|15J] and Mg, Mn, Nb, Zr Jig]. As a general rule for selecting the a dopant, two criteria should 
be taken into consideration: radioactive properties of the nuclide and the incorporation to the host 
lattice of the ion. 

In our particular case, where the detector has to be applied to the search for very rare events, 
supplementary constraints are imposed on the radio-purity of the chemical compound used for 
the doping. Moreover, for cryogenic applications of Te02 crystals, the incorporation of uneven 
ions in the Te02 lattice is to be avoided because of the possible problems of cooling paramagnetic 
materials at very low temperatures. 

The Sm doped crystal reported in the present work was grown by a modified Bridgman method 



described in detail in [17]. For the growth, 6N purity Te02 powder was used. Tellurium oxide raw 



material was synthesized at SICCAS, Shanghai China, in a process described in detail in [18Q. The 
crystal was grown in a dedicated furnace used for R&D purpose where ingots of typically 30 x 25 
x 120 mm can be obtained. The crystal growth followed the protocol described in detail in JI^]. 
Sm dopant was added to the raw material powder in oxide form using Sni203 powder of 6N pu- 
rity. The doping process was performed in two steps. The Sni203/Te02 powder was first melted 
in order to obtain an ingot of homogeneous oxides mixture from which a small sample was taken 
and used as dopant for the final crystal growth with a nominal dopant concentration of 5T0~ 6 g/g 
of Sm203/Te02 . A charge of approximately 400 g (crystal seed excluded) was used for the final 
growth which resulted in an ingot of approximately 65 cm 3 . The grown ingot was colorless and 
free of cracks, bubbles and/or inclusions except for a small region due to a thermal instability dur- 
ing the growth process (see Fig. [[]). It was subjected to special post-growth thermal treatment in 
order to avoid cracks during mechanical processing. A sample of approximately 30x24x28 mm 3 
(m=l 16.65 g) was extracted for cryogenic measurement from the region no.4 of the as grown 
crystal. The sample was cut and X-ray oriented with a precision better than 0.5° following the pro- 



cedure described in [19] except for the chemical etching. Different samples were taken from each 
region ( 1 to 6) for ICP-MS measurements aimed at checking the general quality of the crystal and 
especially the uniform distribution of dopant along the growth direction. ICP-MS measurements 



were made using an 'Agilent Technologies 7500 Series" instrument Q19Q. The measurements were 
made on samarium isotopes free of instrumental interferences: 149 Sm, 152 Sm and 154 Sm. For each 
isotope the Sm concentration was inferred assuming the natural abundance. The mean value is 
found to be of the order of 30 ppb which gives for the 1 16.65 g sample approximately 3.5 /ig Sm. 

3. Experimental set-up and Data Acquisition 

The Sm doped Te02 crystal was operated in a cryogenic setup located deep underground in the 
Hall A of National Laboratory of INFN at Gran Sasso. 



The crystal was equipped with two Neutron Transmutation Doped Ge thermistors (NTD) [ 2C ] 



of 3x3x1 mm 3 , thermally coupled to the crystal surface with 9 epoxy glue spots (~0.8 mm di- 



ameter). The electrical conductivity of NTD, which is due to variable range hopping (VRH) [21] 
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Figure 1. Sm doped TeC>2 crystal as grown. Different samples were cut for analysis. The cloudy (no. 3) 
region is due to a thermal instability during the growth process. 



of the electrons, depends strongly on the temperature. The resistance changes with temperature 
according to [22] 

R = R Q . e \ T ) (3.1) 

and then resistance variation can be used to measure effectively the heat signal produced by inter- 
acting particles. The R and T of the two NTDs were measured to be respectively ~1I2 and ~3K. 
At the working temperature of ~10 mK the value of R is ~ 100 M£l. 

A resistor of ~ 300 kO, realized with an heavily doped meander on a 2.33x2.4x0.6x mm 3 
silicon chip, was glued to the crystal and used as a heater to stabilize the gain of the bolometer. 

The crystal was mounted in an Oxygen Free High Conductivity (OFHC) copper structure and 
kept in position by PTFE tips. The L shape of the PTFE pieces was chosen to profit from the high 
thermal contraction of the Teflon keeping tightly the crystal. 

Finally, the detector was mounted and operated at ~ 10 mK, cooled by an Oxford 1000 
3 He/ 4 He dilution refrigerator. The crystal holder was mechanically decoupled from the cryostat 
in order to minimize noise vibrations induced by the cryogenic facility. A weak thermal coupling 
between the mixing chamber of the dilution refrigerator and the crystal holder was realized by 
means of thin high conductivity copper strips. The holder was also equipped with one NTD ther- 
mometer and one heater. In this way it was possible to stabilize the temperature of the crystal 
holder using a feedback device flffi). 

The cryostat was heavily shielded both internally and externally in order to decrease the y 
background on the detector coming from radioactive materials. Below and above the detector there 



was a ~ 10 cm thick layer of low-activity ancient Roman lead Q24[ |25fl . Around the sides of the 
detector holder an additional 1.2 cm thick cylindrical Roman lead shield was present. The cryostat 
was surrounded by an additional external shielding composed by 10 cm low-activity lead, 10 cm 
standard lead followed by 10 cm borated polyethylene. The latter allows to reduce the neutron flux 
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on the detector thanks to the high efficiency in thermalizing fast neutrons and to the high neutron 
capture cross section for thermal neutrons of 10 B. The entire setup was enclosed in a Faraday cage 
to reduce electromagnetic interference. 

The thermistors were biased through two room temperature 27 G£2 load resistors. The large 
ratio between their resistances and those of the thermistors allows to have negligible parallel noise. 

The read-out of the thermistors was performed through a room temperature DC coupled dif- 
ferential front-end [26] followed by a second stage of amplification both located on the top of the 
cryostat. After the second stage, and close to the DAQ (a 18 bit NI-PXI 6284 ADC unit), a 6 pole 
roll-off active Bessel filter acted as antialiasing filter. 

The entire waveform of each triggered pulse was sampled with a rate of 1 kHz and recorded. 
The typical bandwidth is approximately 10 Hz, with signal rise and decay times of order of 40 and 
200 ms, respectively. 

The two thermistors are read by independent electronic channels and will be denoted in the 
following as Channel 1 and Channel 2. 

A more detailed description of the electronics and the cryogenic facility can be found in 

Ref. 



4. Data Analysis and Results 



Two separate sets of data have been collected at different holder temperatures (runl: ~10 mK, run 
2: ~15 mK ) in order to investigate possible temperature dependencies of the Quenching Factor 
and differences of the pulse shape between a and jS/y particles. 

The energy calibration is performed using four 232 Th r sources inserted inside the external 
lead cryostat shield. The pulse amplitude (A) is estimated by means of an Optimum Filter tech- 
nique [27]. The Channel 1 calibration spectrum in run 1 is shown in Fig. |2[ Gamma lines from the 
232 Th decay chain are visible in the spectrum and listed in Table [|. The 147 Sm a line is also clearly 
visible as well as another a line at 5.407 MeV from 210 Po. This contamination originates during 
the crystal growth. The line from 210 Po is not considered in our QF calculation since it lies in an 
energy region where the calibration is not reliable at the desired level of accuracy. 

The calibration peaks are fitted using a Gaussian function + linear background. 



Isotope - source 


Energy [keV] 


e + e~ annihilation 


511.0 


208 -pi 


583.2 


228 Ac 


911.2 


228 Ac 


968.9 


2614.5 keV DE 


1592.5 


2614.5 keV SE 


2103.5 


208 j j 


2614.5 



Table 1. Calibration peaks from the 232 Th source. 
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5407 



Energy [keV] 

Figure 2. Calibration spectrum of Channel 1 - run 1 . The calibration is performed with external 232 Th y 
sources. The main peaks used for calibration are labeled together with the 147 Sm a line at 2310 keV and the 
210 Po a line at 5407 keV. 



The calibration function is a third order polynomial with zero intercept. The Quenching Factor 
is extracted from a simultaneous fit to the y's and the a line adding the QF as a free additional 
parameter. To properly take into account uncertainties on fitted amplitudes (A,) and transition 
energies (E ( ), the following likelihood function has been minimized: 

1 ( Ei-nAj) \ 2 1 ( Eg~P(A a )/QF \ 2 

-log^ = -Vlog e v I -log-= — e \ ^«« > (4.1) 

where the index runs over the calibration peaks, P(A) is the polynomial calibration function to be 
estimated and a, a is computed as: 



CT '.« = y4« + (li CT ^) (4 - 2) 

The QF calculated for both channels and runs is reported in Table |[ The correlation be- 
tween the QFs has been evaluated on the energy distribution of coincidence events in the a peak, 
recorded simultaneously by the two channels. Taking into account the full covariance matrix, a 
QF = 1.00755 ± 0.00066 and QF = 1.00765 ±0.00066 is found for runl and run2 respectively. 
The error is dominated by the uncertainty on the 147 Sm Q-value. No temperature dependence is 
observed; the average QF is 1.0076 ±0.0005. 

4.1 Energy Resolution and Pulse Shape 

We have performed a study on the dependence of the energy resolution of the y ( 232 Th) and the a 
( 147 Sm and 210 Po) peaks as a function of energy. 



The Optimum Filter [27], used for amplitude estimation, reduces the noise superimposed to 



the signal, maximizing the signal to noise ratio. The ultimate resolution is governed by the noise 
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Run 


Channel 


QF 


QF error 


1 


1 


1.0076 


0.0007 


1 


2 


1.0075 


0.0007 


2 


1 


1.0078 


0.0007 


2 


2 


1.0075 


0.0007 



Table 2. QF calculated for the two channels for both runs. The error is dominated by the uncertainty on the 
147 Sm Q-value (see Introduction). 



at the filter output, and does not depend on the pulse amplitude. We estimated this resolution, Oof, 
on data samples recorded randomly and without triggered events. 

The energy resolution of a and heater peaks has been calculated with Gaussian + linear back- 
ground fits (see for example the fit on the 147 Sm line in Fig. [|). The quantity a 2 is shown as a 
function of energy in Fig. |]. The squared heater resolution o# (black triangle) is well consistent 
with Gq F (dashed line), indicating that no physical process other than the noise contributes to it. 
On the other end, the resolution of a lines lies above Oq F and exhibits a dependence on the energy. 



« 22 
^ 20 
1 18 
§ 16 
14 
12 
10 
8 
6 
4 
2 

2%10 2315 2320 2325 2330 2335 2340 

Energy [keV] 

Figure 3. Gaussian + linear background fit of the 147 Sm line for Channel 1 - run 2. The energy (not corrected 
for the Quenching Factor) is (2328.4 ± 0. 14) keV, with resolution of (0.75 ± 0. 12) keV. 

A linear fit on the a points of Fig. ^ is performed, assuming that the resolution at zero energy 
is equal to the heater resolution (c#) with the corresponding error. The energy dependence of the 
a resolution is well in agreement with the behavior of a calorimeter, where not all the energy is 
eventually collected: 

a [keV] = (0.56 ± 0.02) (0.010 ± 0.002) ^EpseV]. 

This behaviour indicates that in a macrobolometer, where all the energy is expected to be trans- 
formed into heat and measured, some phenomena are responsible for a missing part. We performed 
the same analysis on Channel 2 and found a consistent energy dependence: (0.012 ±0.003)\/keV. 

A more intriguing feature emerges from the energy resolution measured on the 208 T1 y line. 
Its value is worse than the one of an a line of equivalent energy: 

of, = (5.9±2.3)keV 2 

<2614.5keV) = (0.6 ± 0. 1) keV 2 
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0.8 
0.6 
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0.2 


Figure 4. Squared energy resolution (a 2 ) for a's (dots) and heater (triangle) as a function of energy for 
Channel 1 - run 2. The theoretical Optimum Filter resolution {Oq F ) is shown for comparison (dashed line). 
The linear fit on a points is performed assuming that the resolution at zero energy is equal to the heater 
resolution (ojj) with the corresponding error. 

The crystal is however too small and the statistics too limited to allow a detailed study of the energy 
resolutions of the /3/y branch. 

We investigated the difference in pulse shape between a and /3/y particles. In bolometric 
detectors the shape of the signal is not constant with energy as detailed in [Bj. Shape parame- 
ters, therefore, have been compared for events in the 147 Sm a peak ([2320-2334] keV) and in the 
sideband regions ([2200-2320] keV, [2340-2460] keV), mainly populated by /3/y's. 

The rise and decay times, computed as the time difference between the 10% and the 90% of 
the leading edge and the 90% and 30% of the trailing edge respectively, are reported in Table ||. 
The error is dominated by the resolution of the sampling period. No difference is found within the 
error. 



Particle 


Rise Time [ms] 


Decay Time [ms] 


a 


49.0 ± 0.4 


219.6 ± 0.4 


p/y 


49.0 ± 0.4 


219.9 ± 0.4 



Table 3. The rise and decay times for a and /3/y particles in the 147 Sm a region and in the sideband regions. 




5. Conclusions 

We characterized the response function of a Te02 bolometer to a particles with energy close to 
the Q- value of the 130 Te Neutrinoless Double Beta Decay. To achieve this goal we dissolved into a 
Te02 crystal a small amount of nat Sm at the growth phase. 

A Quenching Factor for a particles of (1.0076 ±0.0005) is found, independent of the temper- 
ature in the investigated range. 

The energy resolution on a peaks shows a standard calorimeter energy dependence: a [keV] = 
(0.56 ±0.02) (0.010 ± 0.002) yEjkeVj. Resolution on a ypeak is found to be worse compared 
to an a line of equivalent energy, but the crystal is too small and the statistics too limited to draw 
any conclusion. Signal pulse shape shows no difference between a and /3/y particles. 
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Finally, the presence of a long-lived monochromatic a line in the proximity of the Te Dou- 
ble Beta Decay Rol, could be used to correct thermal gain variation thus allowing less frequent 
calibrations. 
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